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Abstract. We study the late-time evolution of the Universe where dark energy (DE) is 
parametrised by a modified generalised Chaplygin gas (mCCG) on top of cold dark matter 
(CDM). We also take into account the radiation content of the Universe. In this context, the 
late stage of the evolution of the universe refers to the epoch where CDM is already clus¬ 
tered into inhomogeneously distributed discrete structures (galaxies, groups and clusters of 
galaxies). Under these conditions, the mechanical approach is an adequate tool to study the 
Universe deep inside the cell of uniformity. To be more accurate, we study scalar perturba¬ 
tions of the Friedmann-Lemaitre-Robertson-Walker metric due to inhomogeneities of CDM 
as well as fluctuations of radiation and mCCG, the later driving the late-time acceleration of 
the universe. Our analysis applies as well to the case where mGGG plays the role of DM and 
DE. We select the sets of parameters of the mGGG that are compatible with the mechanical 
approach. These sets define prospective inGCG models. By comparing the selected sets of 
models with some of the latest observational data results, we conclude that the mGGG is in 
tight agreement with those observations particularly for a mGCG playing the role of DE and 
DM. 
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1 Introduction 


Modern observational data indicates that the universe started recently to accelerate. 
However, a satisfactory explanation for this acceleration has not been found so far. According 
to observations, something named DE constitutes about 69% of the total value of the energy 
density in the Universe [1]. The nature of this constituent is a challenge for modern cosmology 
and theoretical physics. A number of scenarios has been proposed to explain the problem 
of DE, which can be divided into two groups. On the one hand, the recent speed up of the 
Universe could be explained through consistent modifications of gravity [2-5]. For example, 
the nonlinear f{R) theories are bright representatives of such modified theories (c.f., for 
example. Ref. [ 6 ] for a very recent account on those models). On the other hand, the Universe 
may be filled with a substance which causes the accelerated expansion [7-11]. For example, 
this substance could be a quintessence field or a barotropic perfect fluid. The Chaplygin gas 
is one of the most popular representatives of those fluids [12-14]. 

Indeed, the Chaplygin gas and its generalised version could play the roles of CDM and 
DE [12-14]. Despite this appealing fact, it has been equally shown that the generalised 
Chaplygin gas (GCG), unless a stringent limit is imposed on the parameter of the unified 
model, instabilities and oscillations on the matter power spectrum can arise from the adiabatic 
pressure perturbations of the GCG (c.f. for example [15]). Those instabilities can be avoided 
by considering the GCG unified model as a vacuum-CDM interacting scenario where the 
perturbations are non-adiabatic, being what was named the “geodesic” model in [15] a clear 
representative. Other approaches to overcome the GGG instabilities as a unified model can 
be found in [16, 17]. On the other hand, the GCG could simply play the role of DE and even 
avoid the Big Rip singularity as was shown in Ref. [18]. Furthermore, the GGG can describe 
some primitive epochs of the Universe [19-21] and alleviate the observed low quadruple of 
the GMB [22]. In this work, we will consider the most general type of Chaplygin gas, named 
the modified generalised Chaplygin gas (mGCG) [23], as the substract causing the recent 
speed up of the Universe and in principle unrelated to the CDM component in the Universe. 

Obviously, to explain the late time acceleration of the Universe, the Chaplygin gas 
model or any of its different versions should be valid for a certain period of time, starting 
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from some early time, roughly since red-shift 2 : ~ 0(10), and continuing for some futnre 
period with respect to the present time. In our present work, we consider the case when 
the mGCG induces no future singularities (related to a divergence of the energy density or 
pressnre of the DE component) in the Universe [24-26]. Therefore, the mGGG is valid for 
arbitrarily big values of the scale factor, a. This means that we consider models where the 
scale factor of the Universe does not reach a finite maximum allowed value; e.g. like in a 
Big Freeze singularity [25] We consider the Universe at the late stage of its evolution 
when the galaxies and the group of galaxies are already formed. At this stage, and inside 
the cell of uniformity, the Universe is highly inhomogeneous and the mechanical approach 
[27-29] is an adequate approximation for its description. The theory of scalar pertnrbations 
within the mechanical approach was proposed in [27, 29] . One of the main advantages of this 
approach consists in the fact that it gives the possibility to investigate analytically different 
cosmological models from the point of view of their compatibility with the theory of scalar 
perturbations (see e.g. [30-32]). 

The main idea of such investigations is the following. As we noted above, the considered 
cosmological models are valid for a certain period of time. It is clear that the scalar per¬ 
turbation equations must be consistent during all that period. Goncerning the mechanical 
approach, its accuracy becomes better and better with the growth of the scale factor a. If the 
scale factor of the models is not limited from above, then we can investigate the equations of 
the scalar perturbations for arbitrarily big values of a. This allows us to simplify considerably 
our analysis and look for a set of parameters (if such set exists) compatible with the scalar 
perturbation equations. The models with these sets of parameters are expected to be viable 
as they are consistent with the theory of scalar perturbations for big values of a and at the 
same time can provide a period of acceleration at present. Then, the models selected as 
viable should be tested (e.g. for different ranges of a or redshift) against other observational 
data. If the parameters do not belong to these sets, the corresponding models are nnsuitable 
models in the future. This means that such parametrizations can be used only for a limited 
period of time, starting e.g. from the time of the last scattering surface till present. In our 
paper we find a number of combinations of the parameters for the mGGG models which are 
consistent with the equations of the scalar perturbations in the limit of big scale factors (i.e. 
in the late Universe). These sets define prospective Ghaplygin gas models. 

The paper is structured as follows. In Sec. 2 we review the mGGG cosmology. In Sec. 3, 
we consider the theory of the scalar perturbation for the mGGG within the mechanical 
approach and find the sets of parameters for which the equations of the scalar perturbations 
are consistent. The main results are summarised and discussed in the concluding Sec. 4. 


2 Reviewing the mGCG 


^In the mechanical approach, the functions describing, e.g., the energy density and the pressure for the 
different matter sources as well as their fluctuations are expanded in powers of the inverse of the scale factor, 
1/a. This procedure is mathematically well defined if the scale factor can go to infinity; i.e. 1/a = 0 is the 
point of expansion. On that case, we know for sure that the terms proportional to (1/a)" can be disregarded 
asymptotically as compared with (l/a)"*, if m < n. If the scale factor of the Universe has a maximum size 
in its future expansion, the validity of the mechanical procedure is no longer guaranteed and a separate and 
additional investigation is required. See also footnote 4 and the paragraph related to that footnote. 
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The barotropic perfect fluid called the mGCG is characterized by the equation of state 
(EoS) [23] 

PCh = /3ech - (1 +, (2.1) 

^Ch 

where j3, A and a are the constant parameters of the model. The mechanical approach for 
the particular cases a = — 1 as well as ^4 = 0 were considered in our previous papers [30, 33] 
and we exclude them here. 

Within the parametrization (2.1), the case (3 = —1 is an exceptional one, corresponding 
to the EoS of a cosmological constant. Obviously, if /3 7 ^ —1 the equation of state ( 2 . 1 ) is 
equivalent to the following one: 

Pch = P^Ch - — ) ( 2 - 2 ) 

^Ch 

where A = [\ + I5)A. If we ignore this equality, we can consider (2.2) as an independent 
parametrization (from the EoS (2.1)) where the case (5 = —1 is not reduced to the cosmolog¬ 
ical constant. This particular case will be studied in subsection 3.5 below. 

Let us investigate first the case /3 7 ^ — 1 and a 7 ^ —1. Then, the conservation equation 
results in the following dependence of the energy density on the scale factor a: 

1 

1 + a 

(2.3) 


ech = ech,0 




where oq and ech,o are the current values of the scale factor and the energy density of the 
mGGG, and Ag = A/e]^^. Here, we exclude the case Ag = 0, as it corresponds to ^4 = 0, 
and the case Ag = 1 because it reduces to the trivial cosmological constant case. 

We next analyse the behaviour of the mGGG depending on the values of its parameters 
Ag, a, and /?. If we define the quantities a* and ^ as 


— Q/q 


l-Ag 


Ag 


3(l + /3)(l + a) 


and ^ = (1 + /3)(1 + a), 


we can rewrite Eq. (2.3) as 


ech = ech,o l^sl 


sgn (A^) + sgn (1 - Ag) (^^'j 


1 

l-\-a 


(2.4) 


(2.5) 


Here, sgn( 3 :) is the sign function which returns 1 ( — 1) if x is positive (negative) and 0 if 
X = 0. Notice that since 7 ^ 0,1, a* is a finite positive quantity, and that ^ 7 ^ 0 as long as 
a 7 ^ —1 and /3 7 ^ — 1. 

From Eq. (2.5) we find that in order for the energy density to be well defined, the right- 
hand-side (rhs) of such an equation must be positive. This is ensured by setting (1 -|- a)~^ = 
2 n, where n is an integer, or, for a general a, by requiring that the quantity inside the square 
brackets be positive, 


sgn(^ 5 ) -h sgn(l - Ag) ^ > 0. (2.6) 

When we analyse this inequality, we find the following restrictions for the value of the scale 
factor, summarised in Table 1, 
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1 

1 

1 

0 < a < a* 

0<e 

1 
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0 < a < +00 

1 

1 

1 

O* < O < +00 


Table 1. Intervals of the value of the scale factor, a, for which the inequality (2.6) is satisfied. 


e<o 

0<e 

a* < a < +00 

0 < a < a* 

a, « a ^ ~ ( a ) 

a < a* ^ ££hM ~ Ag 

£ch ,0 ' ' ^ 

a < — 1 lim ech(«) = +oo 

a^a* 

a < —1 ^ lim ech(a) = +oo 

a—>-a* 

— 1 < a =!► lim ech(a) = 0 

a^a* 

— 1 < a => lim ech(a) = 0 

Table 2. Behaviour of the mGCG for < 0. 

e<o 

o<e 

0 < o < +00 

0 < a < +00 

a < a* ^ ~ 1^1^ 

^Ch,0 

a«a, ^ ^ 

a, « a ^ ~ ( a ) 

a* ~ 

a/ 1 ^ 

' ech,o ' ' 

1 ^ = |2Ag|^ 

' ech,o ' ' 


Table 3. Behaviour of the mGCG for 0 < Ag < 1. 


(i) For < 0 the energy density of the mGCG is only defined either to the left or right 
of the value a*, depending on the sign of 

(ii) For 0 < Ag < 1 the energy density of the mGCG is defined for all values of the scale 
factor. 

(hi) For As > 1 the energy density of the mGCG is only well defined to the left or right of 
a*, depending on the sign of 

Having obtained the regions where the energy density of the mGCG is well defined, we 
can study its behaviour near the limits of those regions. Notice that when Ag 0 (0,1) and 
1 + a < 0, the mGCG has a singularity at a = o*. For Ag G (0,1) the mGCG interpolates 
smoothly between a cosmological constant and a perfect fluid with a constant parameter 
of EoS, with a* marking the transition point between these two behaviours. Notice that 
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e<o 

o<e 

0 < a < a* 

a* < a < 00 

a « a* ^ « |H,| 1+- (“;)!+“ 

a* « a ^ ~ ( a‘) 

a < — 1 =k lim ech(a) = +00 

Cl — 

a < — 1 =k lim ech(a) = +00 

Cl — 

— !<«=> lim ech(a) = 0 

— l<a^ lim ech(a) = 0 


Table 4. Behaviour of the mGCG for 1 < Ag. 


the effective cosmological constant behaviour of the mGCG happens in the future (past) 
whenever ^ > 0 (^ < 0) (c.f. Eq. (2.5)). The results obtained are summarized in Tables 2, 
3, and 4. We will restrict our analysis to those cases where the mGCG energy density is 
well defined as a —>■ +oo. For these cases, we obtain late time acceleration whenever ^ > 0, 
where, as stated above, the mGCG mimics a cosmological constant; or when ^ < 0 and 
/3 < —1/3, where the mGCG behaves as a perfect fluid with constant parameter of EoS 
j3. Notice that in the latter we obtain phantom behaviour at late time whenever (5 < —1 
(because the asymptotic effective parameter of the EoS is proportional to j3). 

3 Scalar perturbations: suitable and non-suitable cases 


The theory of the scalar perturbations for the late Universe filled with a barotropic 
perfect fluid with an arbitrary EoS px = Pxi^x) was recently considered in our paper [33], 
where Eqs. (2.17) and (2.18) of that paper are the master equations^. In the case of a mGCG, 
these equations read 


- (^Ch + PCh) = ^PCh + oferad2 j (3.1) 

c^a 3 

A /^c% Kc a KC a 

/\ip 3/C(^ 2 ^ H ^ ^'^rad2 i (^*^) 

where k = SttGtv/c^ (c is the speed of light and Gat is the Newton’s gravitational constant), 

K- = —1,0,+1 for open, flat and closed Universes, respectively, and the Laplace operator A 
is defined via the conformal spatial metrics. Here, y? is the comoving gravitational poten¬ 
tial and 5pc is the difference between the real and average comoving rest mass densities^: 
5Pc = Pc — Pc of inhomogeneities in the form of galaxies and group of galaxies which may or 

^ Fluctuations of an additional perfect fluid (which can be responsible for the late-time acceleration of the 
Universe) arise due to inhomogeneities of dust-like matter (e.g., galaxies). These fluctuations also form their 
own inhomogeneities. In the mechanical approach, it is supposed that the velocities of displacement of such 
inhomogeneities is of the order of the peculiar velocities of inhomogeneities of dust-like matter, i.e. they are 
non-relativistic. In some sense, these two types of inhomogeneities are ’’coupled” to each other [34]. As we 
can see from the master equation below, this results in a change of the gravitational potential. The flattening 
of rotation curves can be one of the consequences of such a change [31]. 

® The comoving rest mass density and the scale factor have dimensions of mass and length, respectively 
[27]. It is worth noting also that the combinations ip/(c^a) and Ka^e are dimensionless. 
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may not include CDM, depending on whether the mGCG plays the role of a unified model or 
just the role of DE (see details in [27, 29]). The analytical results of our work do not depend 
on which of the two background models we are considering. However, the comparison to 
the observations that we do on the Gonclusions applies only to the unified model. Both ip 
and 5pc do not depend explicitly on time, fech and dpch describe fluctuations of the mGGG 
energy density and pressure, respectively. We also take into account the presence of radiation 
in the Universe. This means that to derive Eqs. (3.1) and (3.2), we consider terms up to 
0(l/a^) inclusive and neglect those of the order o(l/o^). This is the accuracy of our study. 
We also consider the possibility that the inhomogeneities and the mGGG fluctuations can 
contribute to the fluctuations of radiation. Therefore, we introduce feradi and (ierad2 which 
are the fluctuations of radiation associated with the matter inhomogeneities and the mGGG 
fluctuations, respectively. Of course, experiments measure the total fluctuation of radia¬ 
tion. Nevertheless, theoretically we can consider these contributions separately. In our paper 
[29], we have shown that deradi = —This relation is already taken into account in 
Eqs. (3.1) and (3.2). With respect to the mGGG, the contribution to the total ferad may take 
place directly via non-zero (5erad2 in Eqs. (3.1) and (3.2). There is also the possibility that the 
mGGG contributes to the radiation fluctuations even if (ierad2 = 0. This can be easily seen 
from Eq. (3.2) where the comoving gravitational potential depends on the presence of dech- 
On the other hand, as we have mentioned above, p defines feradi- Given these two possible 
scenarios, we include (5erad2 into our equations. For these two radiative fluctuation contri¬ 
butions, we have the same equations of state: dpradi = (l/3)feradi) ^PvaA 2 = (l/3)5erad2 ■ 
Therefore, from the matter conservation equation we get (5erad2 ~ 1/a^ (see the appendix in 
Ref. [33]). 

From the mGGG EoS (2.1) we get 


(Jpch 


H,(^-l) + (l-H,)/3(^)3« 
T, + (l-H,)(^)3« 


fcch • 


(3.3) 


Then, taking into account Eqs. (2.1) and (2.3), we get from Eq. (3.1) the expression for fech: 


<^ech = - 


oo \ 3?' 


1 

1 + Q 


+ (1 - H,) (^) 

1 H, + (l-^)(^) 


(1-H,)(1 + /3)(^)3« e-ci,o 
H,(e-1) + (1-H,)/3(M)3« c2a 




3? 


ferad2 


(3.4) 


Now, we can investigate the consistency of Eq. (3.2) with dech given in Eq. (3.4) for 
different values of Before tackling this analysis, we want to make the following remark. 
As we mentioned above, the system of Eqs. (3.1) and (3.2) was derived for the late Universe, 
i.e. for rather big values of the scale factor o, starting e.g. from the present moment. In 
these equations we neglect terms 0(1/0“^) (see footnote 2 in [33]). Obviously, the larger the 
value of a, the more accurate our mechanical approach is. That is, it should work perfectly 
in the future Universe. On the other hand, Eq. (3.4) and consequently (3.2) are considerably 
simplified in the approximation of big values of the scale factor, a, as we are neglecting terms 
of the order o(l/o'^). Following this prescription, we first expand the expression (3.4) in 
powers of 1/a, neglecting terms o(l/a^). Then we substitute this expression into Eq. (3.2) 
and analyse the consistency of the equation obtained. The equation is considered consistent 
if its rhs does not depend on time up to terms 0(l/a) inclusive (bearing in mind that the 
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rhs of this equation has been multiplied by a^). It is clear that the result depends on the 
values of the parameters of the model. If we can define the sets of parameters which make 
this equation consistent, then we could claim that such models are possible candidate for 
prospective Chaplygin gas realisations^. As we can see from Eq. (3.4), the powers of the 
ratio (1/a) are mainly defined by the parameter Therefore, it makes sense to investigate 
this equation separately for different values of 


3.1 ^>I 


In this case, for large values of the scale factor, a, Eq. (3.4) can be approximated as 




Y5erad2 + 0(l/«^) • 


(3.5) 


Therefore, Eq. (3.2) reads 


Aip + 3/C(/9 



Kc^a^ 

2 


1 

3(^-1). 


<5£^rad2 ; 


(3.6) 


where we dropped out the terms o(l/o) (bearing in mind that the rhs of (3.2) has been 
multiplied by o^). The left hand side (Ihs) of this equation does not depend on time®. 
Therefore, the rhs should also not depend on time. Taking into account that (5erad2 ~ l/o^j 
we can achieve this in two ways. Either 6 era.d 2 = 0, which implies that fech = 0 (up to the 
adopted accuracy) and then the mGCG does not cluster, which is physically less reasonable, 
or 


4 ^ 1 - 3/3 

3 ^ ““ 3(l + /3) ’ 


(3.7) 


and the rhs of (3.6) is equal to zero. In this particular case, (Ierad2 can be non-zero, which 
results in non-zero values of (5ech) since (5ech ~ <3erad2 (see Eq. (3.5)). That means the mGCG 
is clustered. 

For example, in the particular case /3 = 0 => a = l/3we get 


PCh — —£Ch,0^s 



and 


ech — ech,o 


As -l- (1 — As) 



(3.8) 


This fluid behaves as dust for small values of a: ech ~ 1/a® (if 0 < As < 1); and as a 
cosmological constant for large values of a: ech ~ const, in full accordance with Table 3 

■^We perform this analysis for big values of the scale factor. It Is clear that the candidates found should be 
investigated further for smaller values of a, e.g. at the present moment where the discrete cosmology approach 
is already valid. We need to compare terms (of the order 0(l/a^)) which were left in Eq. (3.4) with the ones 
we dropped (i.e. o(l/a'^)). Of course, we can do it if we know the observable values of the parameters of the 
model. This is a general remark. However, if a perfect fluid such as the Chaplygin gas is responsible for the 
present day acceleration of the Universe, it is clear that terms 0(l/a®) should not play an important role now. 
According to the observations, they must contribute less than radiation, i.e. terms of the order of 0(l/a'^). 

® We would like to remind that the quantities ip and 5pc are the comoving ones [27]. Therefore, within the 
adopted accuracy when both non-relativistic and weak field limits are applied, they do not depend explicitly 
on time [29]. 
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(right column). Another interesting example corresponds to the choice a = 0 j3 = 1/3. 
Here, 


PCh 


gph — 4Asgch,0 
3 


and 


gCh — gCh,0 


+ (1 ~ ^s) 



(3.9) 


This type of model behaves as a mixture of radiation and a cosmological constant and it has 
been applied to the early Universe in our works [19, 21]. 


3.2 e = l,/3^0 


The case ^ = 1, /3 = 0, implies a = 0 and describes the mixture of the cosmological 
constant and dust, i.e. the ACDM model. Both of these terms (the cosmological constant and 
dust in the form of inhomogeneities) have already been taken into account in the derivation 
of Eqs. (3.1) and (3.2) (see also [27]). Therefore, we exclude this case. 

If we consider the case /3 0 while neglecting terms of order o(l/a'^) in the approxima¬ 

tion of large a, we get from (3.4) 


5gCh ~ 


/3 


1 


1 + 


+ (1 + /3)(1 — ^s) 

A. 


gCh,0 

c^a 


I- As \ao 


^grad2 • 


(3.10) 


Then, Eq. (3.2) takes the form 

1 (3 


A , 1 + /3 ,i+/3«:a2ech,o (l + /3)^ 

-h SfCip - —Spc = -- 2 --^- 


(1-.4.) Affray 


1 As (aQ\-^ Kc^a^ 


3/3 1 - V a / 


ferad2 + 1 1 3 ^ 


a / 
1 \ K(P‘0?‘ 




ferad2. (3.11) 


Because (5ei.ad2 ~ l/a^ and the comoving gravitational potential ip does not depend on time 
(as well as on a), this equation is consistent only if the following system of equations are 
compatible: 


1 + /3 
/3 

(l + /3)^ 
/3 




1 


3/3 1 - A, 


-3 K(?C? 


^gra,d2 0 . 


(l-A,)A^(:^j 


00^3 KO^ech.O 


ip A 


^ 3/3 j 


1 \ K(?(J? 


^grad2 0 . 


(3.12) 

(3.13) 


It takes place e.g. if 5erad2 = ip = However, this is a non-physical result because 
the presence of inhomogeneities in the form of galaxies and group of galaxies in the Universe 
must result in a non-zero gravitational potential; 0. Therefore, we exclude such solution. 
Then, from (3.12) and (3.13) we get respectively 


5grad2 — “3(1 + /3)(1 — As)A^ (~) “^^Ch,0 

V a / Aa 

(3.14) 

3/3 — 1 Vo/ Aa 

(3.15) 
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It can be easily seen that these equations are compatible if /3 = 0. However, this contradicts 
the condition /3 7^ 0. Therefore, the case ^ = 1, /3 7^ 0 is forbidden within the mechanical 
approach. 

3.3 0<^<1 


Here, for big values of the scale factor, a, the mGCG fluctuation (3.4) is approximated 


as 


1 +/3l-As fao\3C echfl , 1 ^ 

A., W 3(130*'“*^ 


and after substituting it in Eq. (3.2) we get 


(3.16) 


Aif + 3ICip - —6Pc = /i i+- 

^ ^ 2 ^ 1-^4 


4.4“ 


no 4^ Ka^£Ch,o 




+ 1 + 


1 


3(1-0 


K(?a? 


2 

(5erad2 • 




(3.17) 


Similar to the previous cases, we need to find conditions under which the rhs of this equation 
is independent of time. Because —2 < 3^ — 2 < 1, the first term on the rhs in Eq. (3.17) 
cannot be compensated by the second one, which behaves as (5erad2 rsj 1/4 [33]. The only 
way to solve this problem is to put d£rad 2 = 0 and to make the first term in Eq. (3.17) 
independent on time. The latter condition takes place if 3^ — 2 = 0 ^ = 2/3. Then, the 

mGCG fluctuation reads 


<5ech — 3(1 + /3) 


1-4 


4.. 


s i 1 + 0 f ^ 


■4,] 


gQh.O 

4 a 


ip. 


(3.18) 


Therefore, the viable set of parameters is 




a = -- 


1 + 3/3 


3(1+ /3) 

A particularly interesting case is/3 = —1/3 a = 0, where 

^Ch + 


PCh = 


and 


^Ch — ^Ch.O 


4s + (1 — 4sj [ — j 


a 


(3.19) 


(3.20) 


It describes a mixture of a perfect fluid with the parameter of EoS —1/3 (e.g. frustrated 
network of cosmic strings) and a cosmological constant. This perfect fluid was considered in 
our papers [30, 33] and is a specific case of the model we used in [22]. It can be easily seen 
from Eq. (3.18) that for this particular case fech = 2(1 — 4s)oQe(2;j^ Q<y9/(c^a^) which exactly 
coincides with Eq. (2.26) in Ref. [30] after the evident substitution eo = (1 — 4s)ech,o- 


3.4 ^ < 0 


Originally, the Chaplygin gas model was introduced into cosmology to describe the late 
time acceleration of the Universe. However, in the case of negative the mGCG behaves as 


- 9 ~ 















vacuum: p(j^ = —ech ~const for small scale factors a, and as some form of matter with the 
energy density ech ~ for large a (e.g. as a dust, radiation and frustrated network 

of cosmic strings or spatial curvature for /3 = 0, /3 = 1/3 and /3 = —1/3, respectively.). 
Therefore, such models are more suitable for the description of the early inflationary era 
rather than the late acceleration as was done shown in Refs. [19-22, 35]. For the sake of 
completeness, we also investigate this case. 

In the case ^ < 0, it is convenient to introduce a new parameter fx = —^ > 0. Then, 
Eq. (3.4) reads 


= 


fflQ \ 


(l-^) + ^(^) 


1 

1 +Ck 


(l-As)(l + /3) /oo\3(i+/3) 


2 ^ 


1 + 

3As{p + l){jf^-{l-As)/3 


+ 


5erad2 • 


(3.21) 


It makes sense to conduct a further analysis, depending on the sign of the parameter /?. 


i) P > 0 


Obviously, for /? > 0 we have 

(^ech ~ —■^5erad2 + o{l/a^). (3.22) 

Therefore, Eq. (3.2) is consistent for /3 = 1/3 when its rhs is equal to zero within the 
accuracy of our approach. Therefore, non-zero values of 5erad2 and (5ech are allowed in this 
particular case. Obviously, there is also the less physically motivated case /3 7^ 1/3 where 
5erad2 = = 0. In this case, the condition fech = 0 means that the Chaplygin gas is not 

clustered. 


ii/ /3 = 0 


For this value of /3, we can write a in terms of /i as a = — (1 + /r), while the mGCG 
fluctuations in Eq. (3.21) are 


<^ech = 


1 - (1 — ^s) ''/ao\3-3M ech,o (1 — ^s) /«o'\3 ech.o 


+- 


As 1 + /i 

1 1 - As /ao\-3^J. 


2 


//(I + p) \ a / c^a 


-if 


1 


1 + /r 


3(1 -|- /i) A 
tent 

4 + 3/r Kc^a? 


ia ) '^^’'"'^2 + 3114 


3(1 -I- p) 


<5erad2 


Then, Eq. (3.2) is consistent if the following system of equations is compatible 

1 


SSrad 2 — 




/ao' 

>3 Ka'^Echfi 

— 

1 o ^ 

V a 

2 


= 0 , 


(3.23) 


(3.24) 


1 1-A, 

I + p As 




3' 

\ a J 




-Se 


rad2 


+ (1 — ^s) 


i/ao\3-3M KO^Ech. 


-ip 


= 0.(3.25) 


As we have already mentioned in subsection 3.2, the solution 5erad2 = (/? = 0 is a non-physical 
one and should be excluded. Then, from equations (3.24) and (3.25) we get respectively 


<5erad2 = 


4 -|- 3p 


(1 -A,)-Vm .ao\3 p _ 

-— 




— ^ech,o 

V a / c^a 


(3.26) 
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and 

<5erad2 = —3(1 — ^s) (~) “^^Ch.O ■ (3.27) 

Vo/ c^o 

The compatibility of these two equations results in the constraint /r(4 + 3/r) = —1. Obviously, 
this constraint is not fulfilled for // > 0. Therefore, the case ^ < 0, /3 = 0 is forbidden within 
the mechanical approach. 

in) /3 < 0, /3 7^ —1 


In this case the Chaplygin gas fluctuations have the following structure 


l + /3,^ , 1 


/3 


(1 - 




n=l 


fech ~ - 

where k = [—(3/fj] is the smallest integer value of the ratio —/3/// and 


ao\3(i+/3) 1 

j 2 ^ Qo"^i'ad2 

a / c^o 6 p 


Ni = - 


71. 


1 - A. 


1 + /? 1 + // 

p 13 


(3.28) 


(3.29) 


First, we consider the condition (5erad2 7^ 0. Because /3 < 0, the (5erad2 terms in (3.28) 
and in (3.2) cannot cancel each other. Therefore, the first expression on the rhs of Eq. (3.28) 
should have the structure const, x {1/a^) + o(l/a'^). A simple analysis shows that this does 
not take place. 

If (5erad2 = 0, then, for consistency of Eq. (3.2), the first expression on the rhs of 
Eq. (3.28) should have the structure const, x (1/a^) + 0(1/0“^). This takes place for (3 = —1/3 
and // > 1/3 implying ^ < —1/3. 


3.5 The Little Sibling case: f3 = —1, a 7 ^ —1 


As we mentioned above, the case /3 = — 1 is a special one. We begin with Eq. (2.2) 

VCh = -ech - , (3.30) 

^Ch 

where A 7^ 0. This EoS generalises our recent work [36] (see also Ref. [37]) and it leads 
to what we named there as the Little Sibling of the Big Rip singularity. Integrating the 
conservation equation for the mGCG we obtain the background energy density 

1 

l + a 

, (3.31) 


ech — ech,o 


\ 3(1+Q!)As 

1 + In ( — ) 

00/ 


where A. = A/e]/^Q. Eor a = 0, we recover the case we studied in Ref. [36]. If we require 
that the energy density is well defined, the rhs of the last equation must be positive. This is 
ensured by setting (1 + a)“^ = 2n, where n is an integer, or, for a general a, by requiring 
that the quantity inside the square brackets be positive If we define the quantity 


a* = oo exp 


1 

3(1 + oi)Ag 


®We will disregard the fine-tuned case (1 -I- a) ^ = 2n. 


(3.32) 
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(1 + (y.'jAs <C 0 

0 <C (1 CX^Ag 

0 < a < a* 

a* < a < +00 

a < —1 ^ lim ech(a) = +00 

a < — 1 lim ech(a) = +00 

a — 

a — 

— 1 < a ^ lim ech(a) = 0 

— l<a^ lim ech(a) = 0 


a^a* 


Table 5. Behaviour of the mGCG for /3 = —1. 


we can rewrite Eq. (3.31) as 


£ch = £Ch,0 


ln( ^ 


3(l+a)A, 


1 

1 + a 


(3.33) 


Therefore, the quantity a* defines the minimum or maximum (depending on the sign of 
(1 + a)As) allowed value of the scale factor for which the energy density of the mGCG is well 
defined. The results obtained are summarized in Table 5. In the mechanical approach, we 
require that the background quantities be well defined for arbitrarily large values of the scale 
factor (see footnote 1 above), therefore we must restrict our analysis to the case (l+a)^^ > 0. 

From Eqs. (3.30) and (3.31), we can easily obtain. 


^Ch +PCh — 

/ 


5pch = 


-1 + 


1 + In 


aAs 


Q \ 3(l+a)As 

aoj 


1 + a 


\ 


V 


\ 3(1+Q:)As 

i + i-'il;) 


<^ech ■ 


(3.34) 

(3.35) 


Then, substituting the previous expressions in Eq. (3.1), we get 
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<5ech = -- 



/\ 3(l+o)^s 

1 + a 

A. 

/ + "'(«) J 

ech,o 


l-aAs + In j 


\3(i+o)^s +a 


1 

o 


1 + In -2^ 

' ao 


3(l+a)As 


3 ~ / \ 3(l+a)As 

1 - aAs + In + 


-6e 


rad2 


ao 


(3.36) 


which for big values of the scale factor, a, the same quantity can be approximated as follows: 


(5e, 


Ch 


-A, 


\ 3(1+Q:)yis 

1 + In ( — ) 

ao/ 


' l + a _ 

ech,o , 1 

9 T 0£rad2 
c^a 3 


(3.37) 


The radiation terms in Eqs. (3.37) and in (3.2) result in a non zero combination (4/3)(5erad2 ~ 
1/a^. Obviously, such combination cannot cancel the first term in (3.37) both for a ^ 0 due 
to the logarithmic behaviour (as the logarithmic functions and the power law functions are 
linearly independent) and, in the case a = 0, because the first term behaves as 1/a. If we put 
<^erad2 = 0, Eq. (3.2) is still inconsistent due to the first term in the rhs of Eq. (3.37). This 
conclusion follows from the simple observation that in this case the term a^Sech depends on 
time both for a ^ 0 and for a = 0. Therefore, the case /d = —1, a ^ —1 is forbidden. 
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^ — (1 -b a)(l + /3) 

5^rad2 

fcch 

late time 

e > 1, e / 4/3 

0 

0 

cosm. constant 

CO 

II 

can be non-zero 

(5eCh ~ <5erad2 

cosm. constant 

C = 2/3 

0 

non-zero 

cosm. constant 

C < 0, /3 > 0, /3 / 1/3 

0 

0 

deceleration 

e < 0, /3 = 1/3 

can be non-zero 

Jech ~ <5erad2 

deceleration 

e <-1/3,/3 =-1/3 

0 

non-zero 

zero-acceleration 


Table 6. Values of the parameters of the mGCG model that are compatible with the mechanical 
approach. 

4 Conclusion and discussion 


In our paper, we have studied the Universe at its late stage of evolution and at scales 
much smaller than the cell of uniformity size which is approximately 190 Mpc [29]. At 
such scales, the Universe is highly inhomogeneous and the galaxies and groups of galaxies 
play the role of the inhomogeneities. We have considered the Universe filled with a modi¬ 
fied generalized Chaplygin gas. Its main properties have been reviewed in section 2. The 
qualitative behaviour of the Universe as well as the asymptotic behavior of the mGCG are 
summarised in Tables 1-5. Additionally, we also took into account the presence of radiation 
in our (perturbative) analysis. 

The main aim of this paper is to investigate the compatibility of the mGGG models 
with the theory of scalar perturbations within the mechanical approach [27, 29]. This ap¬ 
proach is adequate for the late Universe and works well for big values of the scale factor. 
Therefore, we consider the equations of the scalar perturbations in this limit and found the 
sets of parameters which make these equations consistent. As a result, we select the sets of 
parameters of the mGGG models compatible with the mechanical approach. We summarise 
all acceptable cases in Table 6. Depending on the sign of the parameter these cases can be 
split into two groups. Positive ^ provides late time acceleration of the Universe and negative 
^ results in deceleration or zero acceleration of the late Universe. The former group is for us 
the most interesting case because the Ghaplygin gas models have been introduced precisely 
to describe the late acceleration of our Universe. Here, we found three possible choices of the 
parameter ^ which are compatible with the mechanical approach. The cases ^ > 1 (except 
^ = 4/3) and ^ = 4/3 look very unnatural. In the first case the mGCG is not clustered in 
the presence of the inhomogeneities in the form of galaxies, even though so far we have no 
observational evidence that DE cluster. In the second case, the mGCG can be clustered. 
However, these fluctuations do not affect the gravitational potential because the rhs of (3.6) 
is equal to zero. The third choice ^ = 2/3 looks the most promising. Here, the mGCG is 
clustered and affect the gravitational potential. 
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There are a number of papers where the latest observational data were used to get 
constraints on the mGCG models, considered there as a unified CDM and DE model. There 
is also a different approach where mGGG models are applied exclusively to describe DE. 
On the background level the difference between these models consists in the cosmological 
parameter Vt^n- If the mGCG model describes both CDM and DE, then ~ 0.046 and this 
parameter corresponds to the baryonic matter. Such type of models are usually referred to as 
the Unified Dark Energy - Dark Matter scenarios.^ If the mGCG model is responsible only 
for DE, then ~ 0.3 and it describes the baryonic and DM sectors. Here, non-interacting 
DM is inhomogeneously clustered according to the standard picture. In our paper, we did not 
use the parameter for our analysis. Hence, the discrepancy between these two approaches 
does not affect the results in Table 6. However, the experimental restrictions on the parameter 
^ can be different for different approaches. Therefore, we should keep in mind this possible 
discrepancy in ^ based on the observational data obtained for different approaches. 

The former approach where the mGCG models describe both CDM and DE was inves¬ 
tigated, for example, in the papers [43, 44]. The authors of the paper [43] have used the 
Union2 dataset of type supernovae la, the observational Hubble data, the cluster X-ray gas 
mass fraction, the baryon acoustic oscillation, and the cosmic microwave background data 
and found in the Icr estimation that ^ G [0.8492,1.4588]. They considered a flat Universe 
model. More wide la estimation was obtained in [44]: ^ € [0.795,2.032]. Here, the authors 
have used the data for Type la supernovae, baryon acoustic oscillations and for the Hubble 
parameter. The Universe was not supposed a priori to be spatially flat. As we can see, the 
values of the parameter ^ > 1 belong to these intervals. Therefore, both cases with ^ > 1 
indicated in Table 6 are in agreement with these restrictions. However, the most reasonable 
value of ^ from a physical point of view is ^ = 2/3 ~ 0.6667, which is out of these intervals. 
Hence, we can conclude that the status/viability of the mGCG model within the mechanical 
approach is unclear. There are four possibilities. The first, and the most extremal one, is 
that the mGCG model is not viable for the description of DE and CDM. A second possibility 
is that the mGCG model is viable only for a rather limited period of time and cannot be 
applied in the future. That is, the mGCG equation of state (2.1) is not a fundamental one. 
Third, the mGCG is not clustered, similarly to what happens with the cosmological con¬ 
stant. Erom a physical point of view, even though so far we have no observational evidence 
that DE clusters, this is not very natural. However, if it takes place, the mGCG model can 
be viable. There is also a fourth possibility. As we mentioned above, in Refs. [43, 44], the 
mGCG describes both DE and CDM while dust like matter is considered only in a baryonic 
form with the cosmological parameter Dm ~ 0.046. However, we can investigate a different 
scenario where the mGCG is related mainly or exclusively with DE. Here, the pressureless 
matter is a combination of baryonic matter and CDM and Dm ~ 0.3. It is quite possible that 
such choice of Dm changes the la interval for ^ in such a way that the physically reasonable 
value 2/3 will belong to it. Then, such model can become viable. It is also worthwhile to 
investigate the effects that the latest data for the CMB has on the la interval for ^ (similar 

^It is worth noting that the GCG models may suffer from the problem associated with unphysical oscil¬ 
lations or blow-up in the matter power spectrum [38]. This problem in the GCG models can be avoided 
with the help of an unique decomposition into dark energy and dark matter [39]. In addition, the growth 
of inhomogeneities in the CG and GCG models was investigated within the Zel’dovich approximation in the 
papers [40, 41], while the growth of the matter perturbations in the mGCG model was studied at the linear 
level in [42] . It was shown in these articles that the large-scale inhomogeneities evolve in good agreement with 
observations. It is also of interest to investigate a possibility of decomposing the mGCG into DE and DM in 
a similar fashion to the one carried out in [39]. 
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to what was done in [45] for the GCG model). A future investigation will show which of 
these four above mentioned possibilities can be realized. 
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